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ABSTRACT 



Ring diagram analysis of solar oscillation power spec- 
tra obtained from MDI data is performed to study 
the velocity fields within the solar convection zone. 
The three dimensional power spectra are fitted to 
a model with a Lorentzian profile in frequency and 
includes the advection of the wave front by horizon- 
tal flows to obtain the two horizontal components of 
flows as a function of the horizontal wave number 
and radial order of the oscillation modes. This in- 
formation is then inverted using the OLA and RLS 
techniques to infer the variation in flow velocity with 
depth. The resulting velocity flelds yield the mean 
rotation velocity at different latitudes which agrees 
reasonably with hclioscismic estimates. The zonal 
flow inferred in the outermost layers also appears to 
be in agreement with other measurements. A merid- 
ional flow from equator polewards is found to have 
an amplitude of about 25 m/s near the surface and 
the amplitude appears to increase with depth. 



Key words: Sun: oscillations; Sun: rotation; Sun: 
interior. 



(Choudhuri et al. 1995; BrummeU ct al. 1998). 

The high degree modes {I > 150) which are trapped 
in the solar envelope have lifetimes which are much 
smaller than the sound travel time around the Sun 
and hence the characteristics of these modes are 
mainly determined by localized average conditions 
rather than average over entire spherical shell. These 
modes can be employed to study large scale flows in- 
side the Sun, using the time-distance helioseismology 
(Duvall et al. 1993, Giles et al. 1997), Ring diagrams 
(Hill 1988; Patron et al. 1997) and other techniques. 

Ring diagram analysis is based on the study of three- 
dimensional power spectra of solar p-modes on a part 
of the solar surface. If we consider a section of this 3d 
spectrum at fixed temporal frequency, the power is 
concentrated along a series of rings that correspond 
to different values of the radial harmonic number n. 
The frequencies of these modes are also affected by 
horizontal flow flelds suitably averaged over the re- 
gion under consideration. Hence, an accurate mea- 
surement of these frequencies will contain the signa- 
ture of large scale flows and can be used to study 
these flows. Since the high degree modes used in 
these studies are trapped in the outermost layers of 
the Sun, such analysis gives information about the 
conditions in the outer 2-3% of the solar radius. 



1. INTRODUCTION 



The rotation rate in the solar interior has been in- 
ferred using the frequency splittings for p-modes 
(Thompson et al. 1996; Schou et al. 1998). How- 
ever, these splitting coefficients of the global p-modes 
are sensitive only to the North-South axisymmetric 
component of rotation rate. To study the nonaxisym- 
metric component of rotation rate and the meridional 
component of flow, other techniques based on 'local' 
modes are required. Since these velocity components 
are comparatively small in magnitude they have not 
been measured very reliably even at the solar surface. 
Apart from these nearly steady flows there could also 
be cellular flows with very large length scales and life- 
times, i.e., the giant cells, which have been believed 
to exist, although there has been no firm evidence 
for such cells (Snodgrass & Howard 1984; Durney et 
al. 1985; Howard 1996). Study of such large scale 
flows is important for understanding the theories of 
solar dynamo and turbulent compressible convection 



2. THE TECHNIQUE 



In this work we have used data from full-disk 
Dopplergrams obtained by the SOI/MDI instrument 
on board SoHO. Selected regions of Dopplergrams 
mapped with Postel's projection are tracked at a rate 
corresponding to the photospheric rotation rate at 
the center of each region to fllter out the photospheric 
rotation velocity from the flow fields. This allows us 
to study the smaller components of the flow which are 
not very well determined from other studies. For each 
tracked region, the images are detrended by subtract- 
ing the running mean over 21 neighboring images to 
filter the series temporally. The detrended images are 
apodized and Fourier transformed in the two spatial 
coordinates and in time to obtain the 3d power spec- 
tra. We have chosen the spatial extent of the region 
to be about 15° x 15° with 128 x 128 pixels in heho- 
graphic longitude and latitude giving a resolution of 
0.03367 Mm-i or 23.437 i?Q^ Each region is tracked 
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Figure 1. Sample (logarithmic) power spectra as a func- 
tion of kx and ky at a fixed frequency. The top panels are 
for the region centered at equator and Carrington longi- 
tude of 60°, at frequencies of around 3 mHz (left) and 4 
mHz (right). The bottom left panel is the spectrum around 
4 mHz for a region centered at 4:0° N latitude and longi- 
tude of 60° . The bottom right panel is the summed spec- 
trum around 4 mHz for regions centered at the equator. 
The gray scale is marked with logarithm of the power. 

for 4096 minutes giving a frequency resolution of 4.07 
/iHz. To minimize effects of foreshortening all the re- 
gions were centered on the central meridian. The 
spectra have been obtained using the usual tasks in 
MDI data processing pipeline. Fig. 1 shows some of 
the sections at constant v of these spectra. 

To extract the flow velocities and other mode param- 
eters from the 3d power spectra we fit a model of the 
form 

exp(Ao + (fc-fco)Ai+A2(i^)"+A3^) 

{v-ckP-U^k^^-Uykyy-'+w'^ y I 

where k"^ = k^ -\- ky, and the 11 parameters 
Aq, Ai, A2, Ari,c,p,Ux,Uy,w, B2 and are deter- 
mined by fitting the spectra using a maximum likeli- 
hood approach (Anderson et al. 1990). Here ko is the 
central value of k in the fitting interval. The A2 and 
A3 terms account for the variation of power along the 
ring. The minimization has been performed using a 
quasi-Newton method based on the BFGS formula 
for updating the Hessian matrix (Antia 1991). The 
form given by Eq. (1) is slightly different from what 
is used by Patron et al. (1996), because we have as- 
sumed some variation in amplitude along the ring, 
coming from the A2 and terms. These were intro- 
duced because the power does appear to vary along 
the ring and the fits in the absence of these terms 
were not satisfactory. 



We fit each ring separately by using the portion of 
power spectrum extending halfway to the adjoining 
rings. For each fit we use a region extending about 
zblOO/iHz from the chosen central frequency. We 
choose the central frequency for fit in the range of 
2 5 mHz as the power outside this range is not sig- 
nificant. The rings corresponding to < n < 6 have 
been fitted. This gives us typically 800 'modes', all 
of them may not be independent as there is a consid- 
erable overlap between adjacent fitting intervals. 

The fitted and Uy for each mode represents an av- 
erage over the entire region in horizontal extent and 
over the vertical region where the mode is trapped. 

We can invert the fitted Ux (or Uy) to infer the varia- 
tion in horizontal flow velocity (or Uy) with depth. 
We use the Regularized Least Squares (RLS) as well 
as the Optimally Localized Averages (OLA) tech- 
niques for inversion. For the purpose of inversion 
the fltted values of Ux and Uy are interpolated to the 

nearest integral value of k (in units of -Rq^) and then 
the kernels computed from a full solar model with 
corresponding value of degree i are used for inver- 
sion. 

We have selected the regions centered at Carrington 
longitudes of 90°, 60°, 30° for rotation 1909 and at 
360°. 330°, 300° for rotation 1910 corresponding to 
period from about May 24 to June 7, 199G. For each 
longitude we select regions centered at latitudes of 0° , 
±10°, ±20°, ±30°, ±40°, ±50° and ±60°. There is 
some overlap between different regions. Apart from 
these individual spectra we have also analyzed spec- 
tra obtained by summing all 6 spectra for a given 
latitude to study the properties averaged over differ- 
ent longitudes. Because of averaging, these spectra 
have better statistics and the error estimates are also 
lower. These averaged spectra can be expected to 
give the average velocity over the range of longitudes 
considered. Most of the inferences in this work have 
been obtained using these averaged spectra. We flt 
each spectrum to obtain the mode parameters includ- 
ing Ux and Uy. 



3. RESULTS 



Following the procedure outlined above we fit the 
form given by Eq. (1) to a suitable region of a 3d 
spectrum. Although other quantities may also be of 
some interest, in this work we restrict our attention 
to the two horizontal components of velocity (Fig. 2) 
obtained by fitting the spectra. The fltted veloci- 
ties for each 'mode' are then inverted to obtain the 
variation of horizontal velocity with depth. Only the 
region r > 0.97-Rq is sampled by the modes used in 
this study and hence the inversions are restricted to 
this region. Some of the results obtained using the 
RLS and OLA techniques are shown in Fig. 3. 

From the inversion results it appears that the lon- 
gitudinal component (ux) is dominated by the aver- 
age rotational velocity. This is a result of the fact 
that tracking is done at the surface rotation rate 
at the centre of the tracked region and hence does 
not account for the variation of the rotation rate 
with depth. The average variation with latitude, 
9 of these components can be determined from the 
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Figure 2. The fitted velocity for summed power spectra 
at different latitudes. In each panel the crosses mark 
the fitted velocity for northern hemisphere while the open 
squares mark that for the southern hemisphere. In the 
right panels the crosses marking northern latitudes gen- 
erally fall in the upper half of the figure, while the squares 
marking southern latitudes fall in the lower half. The lat- 
itudes are marked in the left panel. Error bars are not 
shown for clarity. 
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Figure 3. A sample set of inversions are shown for lati- 
tudes ±40° and longitude 300°. The circles and triangles 
are the OLA inversions for 40° N and 40° 5 respectively. 
The continuous and dashed lines are the RLS results for 
40° TV and 40° 5 respectively, with the dotted lines denoting 
the la error limits. 
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Figure 4- The average horizontal velocity at different lat- 
itudes are compared with the rotation velocity ( after sub- 
tracting out the surface rotation rate used in tracking each 
region) obtained from inversion of splitting-coefficients 
(continuous line). The results obtained using RLS in- 
versions are shown by short-dashed (northern latitudes) 
and long-dashed lines (southern latitudes) with dotted 
lines marking the la error limits. Similarly, the crosses 
(north) and squares (south) with error bars represent the 
results of OLA inversions. 



summed spectra and a reasonable agreement between 
Ux and the rotation rate determined from splitting 

coefficients of the global p-modcs can be seen from 
Fig. 4. The rotation velocity at each latitude can be 
decomposed into the symmetric part [{un + us)/2] 
and an antisymmetric part [{un — us)/2]. The sym- 
metric part can be compared with the rotation ve- 
locity as inferred from the splittings of global modes 
(Antia et al. 1998) and the results arc shown in Fig. 5. 
Since the inversion results using global modes are not 
particularly reliable in the surface regions, the veloc- 
ity profiles obtained from the ring diagram analysis 
supplement those results and support the earlier con- 
clusions. 

Following Kosovichev & Schou (1997) it is possible 
to decompose the rotation velocity into two compo- 
nents, a smooth part [polynomial in terms of cos(^), 

cos^{9) and cos^{9)] and the remaining part which 
has been identified with zonal flows. The zonal flow 
so estimated is shown in Fig. 6. Despite poor lat- 
itudinal resolution in our results the inferred pat- 
tern near the surface is in reasonable agreement with 
the average zonal flow estimated from the splitting 
coefficients for the f-modes from the 360 day MDI 
data. But at deeper depths the pattern changes sig- 
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Figure 5. The average horizontal velocity at different lat- 
itudes (dashed line for RLS and crosses for OLA), i.e., 
the latitudinally symmetric part [[un -\-Ug) /2], compared 
with the rotation velocity (after subtracting out the sur- 
face rotation rate used in tracking each region) obtained 
from inversion of splitting-coefficients (continuous line). 



nificantly and the errors arc also larger. Hence, it 
is not clear if the zonal flow penetrates below about 
7Mm (O.Oli?©) from the surface. 

The difference between rotation velocity at the same 
latitude in the North and South hemispheres is small 
and thus the antisymmetric component of rotation 
rate is not very significant. Some of the difference 
may also be due to some systematic errors in our 
analysis. For example, due to differences in the angle 
of inclination for the regions at same latitude in the 
two hemisphere, the effect of foreshortening will be 
different in the two hemispheres. The antisymmet- 
ric component of the velocities is shown in Fig. 7. 
In particular, it can be seen that at low latitudes 
where the results are more reliable, the antisymmet- 
ric component is rather small, more or less within 
error estimates. 

The latitudinal component of the velocity appears 
to be dominated by the meridional flow from equa- 
tor polewards. The average latitudinal velocity for 
each latitude is shown in Fig. 8, while Fig. 9 shows 
the same as a function of latitude at a few selected 
depths. There is a significant variation in this veloc- 
ity with depth at high latitudes. Since the measure- 
ments are not particularly reliable at high latitudes 
it is difficult to say much about the general form of 
the flow velocity with latitude and moreover it will 
also depend on depth. However, if we fit to a form 
(cf., Giles et al. 1997) 
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Figure 6. The zonal flow, i.e., the rotation velocity after 
removing a smooth component, at a layer just below the 
solar surface. The dotted line and points are RLS results, 
while the dashed line and points are OLA results. The 
error bars shown in this figure represent the error in ro- 
tation velocity and do not include any contribution from 
the smooth component which is subtracted to obtain these 
values. Hence the errors are likely to be underestimated. 
The continuous line represents the average zonal flow ve- 
locity as inferred from the f-modes using the 360 day MDI 
splitting coefficients. 
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Figure 7. The antisymmetric component [{un — us)/2] 
of the rotation velocity plotted as a function of depth for 
various latitudes. The lines are RLS and crosses are OLA 
results. 



then we find mean flow 02 of about 25-40 m/s de- 
pending on the depth, while ai is small, more or less 
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Figure 8. Meridional velocity at different latitudes plot- 
ted as a function of depth. The results obtained using 
RLS inversions are shown by continuous lines (northern 
latitudes) and dashed lines (southern latitudes) with the 
dotted lines showing the la error limits. Similarly, the 
crosses (north) in the positive half and squares (south) 
in the negative half of each panel represent the results of 
OLA inversions. 
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Figure 9. Meridional velocity at different depths plotted 
as a function of latitude. These results have been obtained 
using the OLA technique for inversion. 



comparable to error estimates. It is not clear if the 
assumed form indeed fits the measured variation as 
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Figure 10. Horizontal flow velocities in horizontal planes 
at various depths obtained after the rotation and merid- 
ional velocities (as shown in Fig. 6 & 8) are subtracted 
from the individual measurements. The panels (a), (b), 
(c) and (d) respectively show the residual flow pattern at 
radial distance of 0.995, 0.990, 0.985 and O.QSOi?© as 
inferred by the OLA inversion. The arrow at the top 
marks the scale. The errors in these measurements are 
not shown but they are typically 2-5 m/s depending on 
latitude and depth. 



the resulting is fairly large. The expected decrease 
in the magnitude of Uy at high latitudes, is not very 
clear from the results, though at layers immediately 
below the surface the meridional velocity does ap- 
pear to decrease with latitude at high latitudes. The 
amplitude is about 25 m/s at a radial distance of 
O.997i?0 which is comparable to the estimated value 
obtained by Giles ct al. (1997) from time-distance 
analysis and by Hathaway ct al. (1996) from direct 
Dopplcr measurement at solar surface. The ampli- 
tude appears to increase rapidly with depth to about 
40 m/s, and it is not clear if this form of latitudi- 
nal variation is actually valid at deeper depths where 
the velocity does not show any turnover at high lati- 
tudes. There is no evidence for any change in sign of 
the meridional velocity up to a depth of O.OSi?© or 
21 Mm that is covered in this study. 

In order to extract other components of large scale 
flow fields and to study their variation with latitude 
and longitude we subtract out the average Ux and 
Uy as determined from the summed spectra, from 
those for individual regions. The residual velocities 
are shown in Fig. 10. Some of this velocity could be 
due to supergranules. It may be noted that the ve- 
locities at different longitudes in this figure gives the 
average values centered at different times. The rms 
residual velocity at each of these depths is found to 
be 12 16 m/s. There docs not appear to be any clear 
pattern in this velocity thus suggesting that the giant 
cells if they exist have velocities less than about 10 
m/s, or their life-time is smaller than the duration of 
about 15 days for which the data has been analyzed 
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in this work, or their longitudinal (latitudinal) size is 
less than about 30° (10°). 



4. CONCLUSIONS 



The average longitudinal velocity agrees reasonably 
well with the rotation rate inferred from inversion of 
global p- modes. The zonal flow velocity in the outer- 
most region also agrees with that estimated from the 
splitting coefficient for the global f-modes, though it 
is not clear if this flow continues in deeper layers. The 
antisymmetric component of the rotation velocity is 
small (< 5m/s) and it is not clear if it is significant. 

The dominant signal in the meridional velocity is the 
meridional flow which varies with latitude and has a 
maximum magnitude of about 40 m/s. The ampli- 
tude of meridional component increases from about 
25 m/s just below the surface to 40 m/s in deeper lay- 
ers. The form of variation with latitude is not very 
clear, but the velocity increases with latitude until 
about 40° latitude at all depths. There is no change 
in sign of meridional velocity with depth up to 21 
Mm. 

The residual after removing the dominant rotation 
and meridional flow component has a magnitude of 
about 10 m/s and may represent flows due to the 
giant cells or some residual contribution from super- 
granules. From the absence of any clear pattern in 
these flows it appears that the giant cells if they exist 
have a velocities less than 10 m/s, or have lifetimes 
smaller than 15 days, or their longitudinal size is less 
than about 30°. 
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